Renormalisation group evolution for the AS = 1 effective Hamiltonian with Nf = 2 + 1 
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We discuss the renormalisation group (RG) evolution for the AS = 1 operators in unquenched 
QCD with Nf = 3 (m u = m<j = m s ) or, more generally, Nf =2 + 1 (m u = m<j m s ) flavors. In 
particular, we focus on the specific problem of how to treat the singularities which show up only 
for Nf = 3 or Nf — 2 + 1 in the original solution of Buras et al. for the RG evolution matrix at 
next-to- leading order. On top of Buras et al.'s original treatment, we use a new method of analytic 
continuation to obtain the correct solution in this case. It is free of singularities and can therefore 
be used in numerical analysis of data sets calculated in lattice QCD. 
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I. INTRODUCTION 

In the standard model, the direct CP violation param- 
eter e'/e and weak decays of hadrons are described by a 
low energy effective Hamiltonian which can be obtained 
by decoupling the heavy quarks and gauge bosons. In 
this paper we consider the AS = 1 effective Hamilto- 
nian [J, [3, H, H § which governs e'/e and the AI = 1/2 
rule. The Hamiltonian is composed of hadronic matrix el- 
ements of four fermion operators with Wilson coefficients 
which arc known up to next-to-leading order (NLO) in 
perturbation theory 0, @ • Since the energy scale in kaon 
decays is about 500 MeV, the hadronic matrix elements 
arc dominated by the strong interaction, QCD. Hence 
we must introduce a non-perturbative tool such as lat- 
tice gauge theory in order to calculate them 0, H, H [H| • 

Once we calculate the hadronic matrix elements on the 
lattice, we need to convert them into the correspond- 
ing quantities defined in a continuum renormalization 
scheme such as naive dimensional regularization (NDR); 
this is often called "matching" . When we match the lat- 
tice results to the continuum, we must introduce a match- 
ing scale q* [ll|. A typical choice of q* lies in the range 
from l/a to ir /a 12]. Once we match the lattice results to 
those in the continuum NDR scheme at the q* scale, we 
have two options to combine the Wilson coefficients with 
the continuum results. One is to run the hadronic matrix 
elements from q* down to m c (Nf = 3) and the other is 
to run the Wilson coefficients from m c up to q* . To do 
this we use the RG evolution equation at NLO, which is 
explained in great detail in 0,11, [![. In this paper we 
focus on the RG evolution equation and its solution at 
NLO. 

For three sea quark flavors (Nf — 3 or Nf = 2 + 1), 
singularities arise in the NLO solution given in [l|, even 
though the full RG evolution matrix is finite. This makes 



'Electronic address: 
t Electronic address: 



dadams@phya.snu.ac.kr 



wlee@phya.snu . ac . kr | 



it impossible to calculate the RG evolution matrix numer- 
ically in this case, which is an essential step for the lat- 
tice evaluation of e'/e and kaon decay amplitudes. Two 
unsatisfactory approaches to dealing with this problem 
have been attempted previously in the literature. In 
Ref. 0,0, the Wilson coefficients for N f = 3 were com- 
bined with the hadronic matrix elements calculated using 
the RG evolution matrix with Nf = (quenched QCD). 
In Ref. [1, HH , the singularities were removed artificially 
by putting in an arbitrary cut-off of « 1000 by hand 
in the calculation of the Nf = 3 RG evolution matrix 
(i.e. singular matrix elements were replaced by finite 
ones with this value). 

In this paper we provide the correct solution for this 
problem. Singularities do not arise, and it can therefore 
be used for numerical calculations. The applicability of 
our result is broad enough that any lattice calculation 
regardless of fermion discretization can take advantage 
of it. In fact, the results of this paper are already being 
used for the ongoing data analysis of the staggered e'/e 
project 0,[T3|. 

Unquenched lattice simulations with Nf = 2 + 1 
sea quark flavors are currently underway with a variety 
of fermion discretizations: AsqTad staggered fermions, 
HYP staggered fermions, Wilson clover fermions, twisted 
mass Wilson fermions, domain wall fermions, and over- 
lap fermions (see, e.g., [HI and references therein). These 
simulations are currently focusing on decay constants (/ w 
and /if), hadron spectrum, the indirect CP violation pa- 
rameter Bk, kaon semileptonic form factors, and kaon 
distribution amplitudes 15]. In the coming future, the 
simulations will be extended to address kaon physics such 
as e'/e and the AI =1/2 rule. The results of this paper 
will be needed in connection with this. 

This paper is organized as follows. In Sec. fill we review 
the RG evolution results originally presented in 0] and 
raise the serious problems for Nf = 3. In Sec. IIII| we 
provide the correct solution to the QCD part of the RG 
evolution equation, in which singularities do not arise. In 
Sec. |lVl we go on to present the correct solution to the 
QED part of the RG evolution equation, in which there 
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are also no singularities. We close with some conclusions. 

II. REVIEW OF RG EVOLUTION 

The AS = 1 effective Hamiltonian for non-leptonic 
decays may be written in general as 



(lF Q T {n)-C{n) (1) 



where the index i runs over a basis for the contributing 
operators; in our example this is the basis Qi, Q 2 , • ■■ , 
Qio of of Buras et al. defined in Section 2 of Ref.pj. The 
renormalization group equation for C(/i) is 



C = 1 T (g,a)C 



(2) 



where (Q)^ and (Q) denote the unrenormalized and 
renormalized Green functions, respectively. Z q is the 
quark field renormalization constant and Z is the renor- 
malization constant matrix of the operators Q. The 
anomalous dimension matrices are defined by 



7(5, a) = Z 1 -r^—Z 
am// 



(9) 



which includes QCD and QED contributions. For the 
case at hand, j(g, a) can be expanded in the following 
way (with a = f^): 



l(g,a)= ls (g 2 ) + —T(g 2 ) + 



(10) 



The QCD part of the anomalous dimension, 7 S , can be 

2 

expanded as (with a s = 4-) 



where /3(g) is the QCD beta function: 

,3 



P{g) = -A) 



9 



16tt 2 



ft 



9 



(16tt 2 ) 5 



le 



e 2 g 3 
(16tt 2 ) 2 



with 



A = 11-3/ 



and / 



ft = 102- f/ 



(3) 



d 
f 4 
(4) 

d denoting the number of active flavors, 
u and d being the number of u-type and d-type flavors 
respectively. According to Buras et al., the the contri- 
bution from the/3i e term is negligible and they dropped 
it in Ref.'s [l|, Q. We will include it here for the sake 
of completeness. The j(g, a) matrix is the full 10 x 10 
anomalous dimension matrix, which is given in [l[. 

The solution of the RG equation in Eq. @ for the 
Wilson coefficient functions is given by 



C{n) = U{n,nw>ot)C{p,w) 



(•5) 



ls[9 ' Air ls + (4^) 2% + 1 j 

The QED part of the anomalous dimension, T, can be 
expanded as 

r(ff 2 ) - 7< 0) + ^7« + ■■■ (12) 

The general RG evolution matrix U(mi, m 2 , a) of Eq. ([6|) 
may then be decomposed as follows 



U(mi,m 2 ,a) — U{m\, m 2 ) + 7- R(mi, 1712) (13) 

47T 

ff(mi) ^ T (n')\ 
(m 2 ) P{9 ) J 

l U(mi,m')T T {g , )U(m',m 2 ) 
R(mi,m 2 ) = / dg —- 

Jg{m 2 ) P\9 ) 

(15) 



U(m 1 ,m 2 ) = T g ex.p(j 



The coefficients at the scale nw — 0(Myy) can be eval- 
uated in perturbation theory. The evolution matrix U 
then includes the renormalization group improved per- 
turbative contributions from the scale fiw down to \i. 
For mi < J7i2, 



U(mi, m 2 , a) = T g exp 



9(mi) , n T {g'^) 



lg(m 2 ) 



dg 



/%') 



(6) 



where g is the QCD coupling. Here, T g denotes ordering 
in the coupling constant such that the couplings increase 
from right to left. Note that for g\ ^ g 2 , 



[7(00,7(02)] 7^0 



(7) 



The evaluation of the amputated Green functions with 
insertion of the operators Q gives the following relation 



(Q) 



(0) 



zr 2 



Z(Q) 



(8) 



Here, U(mi, TII2) represents the pure QCD evolution and 
R(mi, m 2 ) describes the additional evolution in the pres- 
ence of the electromagnetic interaction. The leading or- 
der RG equation describing the QED evolution was first 
discussed in [jq |. 

The QCD evolution matrix which Buras et al. pro- 
vided originally can be expressed up to NLO as 



_ a s (m 2 ) 



47T 



(16) 



where U^°'(mi,m 2 ) denotes the evolution matrix in the 
leading logarithmic approximation and J summarizes the 
next-to-leading correction to this evolution. Additional 
terms proportional to a 2 in U(m,\,m 2 ) which do not 
come from U 1 ^ (m\ 1 m 2 ) should be consistently dropped 
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at NLO. Taking V to be a matrix which diagonalizes 
7s , we define the following: 



(0) = V -l 7 (0)Ty 



lb 
G 



y-l 7 (l)Ty 



(17) 
(18) 



where 7 jj denotes a diagonal matrix whose diagonal el- 
ements are the components of the vector 7*' J . Then, 



C/ (0) (mi,m 2 ) = V 



a s (m 2 ) 
a s (mi) 



V~ 



J D 



with 



. 7 1 0) 
° = Wo 
For the matrix J, we find 

J = VSV' 1 

where the elements of S are given by 

(o) Pi Gj 



5 ijTi 



2/?o 2 2/3 + 7l (0) - 7 



(o) 



(19) 



(20) 



(21) 



(22) 



where is a component of and Gy denotes the 



1811. This is the result obtained 



elements of G in Eq. (j 

by Buras et al. [l|, Q- However, when / = 3(= Nf), 
i = 8 and j = 7, the denominator of Eq. (|22|) vanishes: 



2/3o + 7» ( ° ) ~ 7-"' = 0. Hence the solution for the RG 
evolution matrix at NLO has singularities in this case, 
despite the fact that the full RG evolution matrix must 
be finite. This is one of the main issues we will address 
in this paper. This problem is briefly mentioned in [3, [EJ 
without any solution given. The correct treatment for 
this case, which is given in the next section, eradicates 
the singularities in Eqs. (fT6|) and ([22]) . and leads to a 
finite expression for the RG evolution matrix at NLO. 

Next, we turn to the QED part of the evolution matrix 
R(mi,m 2 ) given in Eq. (| 15[) . We can expand R(mi,m 2 ) 
in powers of g 2 as follows: 

%i,m 2 ) = i? (0) (mi,TO 2 ) + R {1) (mi,m 2 ) + ■ ■ (23) 

where R^ is of the order of g 2% . It turns out to be con- 
venient to introduce the matrix K(mi,m 2 ) to represent 
R(mi,m 2 ) as follows: 



(0) 



R(mi, m 2 ) 
K(mi,m 2 ) 

i? (0) (mi,m 2 ) 
R^(mi, m 2 ) 



-^VK(m ll m 2 )V- 1 
Po 

if (0) (mi,m 2 ) 



1 3 

47T ^— ' 



(1) 



(mi, ma) 



VK^im^m^V- 1 



(24) 

(25) 
(26) 



2tt 
"po 

3 

A, 4^ 



The leading order term can be obtained by straightfor- 
ward integration 



(ifW(mi,m 2 )) i 



ill 



(0) 



a,- — o,- — 1 



a s (m 2 )\ a 



a s (mi) 



a s K) / a s (mi) 



(28) 



5 (m 2 ) 



where the a^'s are the components of a in Eq. (|2U)l . and 
the M(°) matrix is given by 



m<°> = v-H 0)T v 



(29) 



Similar to the case of Sij in Eq. (|22|) , there is a singularity 
in Eq. (HHJ) for the clement (7,8) of (K ( - \m 1 ,m 2 )) since 
127 = as + 1 when / = 3(= -/Vf). However, the expression 
in the numerator also vanishes in this case and so the 
singularity is removable. In this case of — aj + 1, 
direct integration leads to the following formula [lj]: 



JT (0) (mi,m 2 ) 



1 



M (0) 

13 a s (mi) 



a s (m 2 ) 
as (mi) 



x In 



a s (mi) 
a s (m 2 ) 



(30) 



The next leading corrections to the QED part of the 
evolution matrix are represented by (mi, m 2 ). We 
introduce 



and 



-(1) = A (0)T_ (0)T /o^ 

P0 P0 



M« = y-i(r( 1 ) + [ 7 (°) T ,J])V 



(32) 



The matrices (mi, m 2 ) are then given as follows: 
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(^ 1) (m 1 ,m 2 ) 



Qij — 



/ a B (m 2 ) \ aj _ ( a s (m 2 ) \ ' 
\«.(nil) J ^o,(mi) J 

In 



\a.(mi) j \a 



■ M 
s(rn 2 ) 



if i + 3 
if i = j 



(33) 
(34) 



K. 



2^(mi,m 2 ) = — a s (m.2)K^\m.i 1 m 2 )S (35) 
/4 1) (mi,m 2 ) = a s (m 1 )SK^(m u m 2 ) (36) 

As one can see in the above equations, all of the 
(tii , Tn%) matrices include S or J. Since S or J 
is singular for / = 3(= Nf), all the K^' matrices are 
also singular. However, the full RG evolution matrix 
R(mi,m 2 ) must always be finite. The correct treatment 
of R(mi,m2) at NLO in this case, given in Section HVl 
modifies the formulae given in Eq. (|33M36[) such that all 
the singularities are eradicated. 



III. HOW TO HANDLE REMOVABLE 
SINGULARITIES FOR N f = 3 (QCD PART) 

For / = 3 (three dynamical flavors), when i = 8 and 



] = 7, A) = 9, 7, 



(0) 



16 and 7 W = 2. Hence, 2/3 



7 



" 7^ = corresponds to a pole in Sij of Eq. (f2"2")l . 



which does not exists in the full evolution matrix and so 
should not appear in the correct solution at NLO. To find 
the latter, we start by rearranging the NLO expression 
for the QCD RG evolution matrix given in Eq. (fTB)) as 
follows 

U(m 1 ,m 2 ) = Uo{m l7 m 2 ) + -j-VA(m 1 , m 2 )V^\37) 

47T 

where 

VA(mi,m2)V~ 1 = a s (mi)JU (mi,m 2 ) 

-a s (m 2 )U Q (mi,m 2 )J (38) 

When S and J matrices are non-singular, the A matrix 
is readily found from Eqs. ([ig ]) -([2"2" ]l and Eq. (J35D to be 

given by 



A{,j Sij 



a s {mx) 



a s (m 2 ) 
a s {mi) 



a s (m 2 ) 



a s (m 2 ) 
as (mi) 



(39) 



When Sij is singular (i.e. i — 8 and j — 7), this expres- 
sion diverges and therefore cannot be used in numerical 
calculations. In this case, the correct finite expression for 
the A matrix is 

Gij (a s (m 2 )\ a% fa s (m 2 )\ 

A v = T7jra s (m 2 ) — In — '- 40 

2/3 \a s (mi) J \a s (mi) J 



The derivation is as follows. For i =/= j , 

s = 

w 2^,(1 + ^-0^) 



(41) 



We regularize the singularity that occurs when cij = a^+l 
(i = 8 and j = 7) by introducing an e shift of aj so that 



aj = at + 1 + e. Then 



(42) 



and Aij is given by 



A, 



-0(e 2 ) 



1 / a s (m 2 ) 

-a s (m 2 ) — r 

e \a s (rni) 



e In 



u s (m 2 ) 
a s (mi) 



(43) 



In the limit of e = 0, we get the claimed result of Eq. (|40| . 
The expression Eq. (|40ll is finite and can therefore be used 
for numerical calculations. In the appendix we give an 
alternative derivation of this result, similar to the pre- 
vious approach of Buras et al. in the non-singular case 



IV. HOW TO HANDLE REMOVABLE 
SINGULARITIES FOR N f = 3 (QED PART) 

Now let us turn to the QED part of the evolution ma- 
trix for / = 3(= Nf). Basically, we want to perform 
the integration in the right-hand side of Eq. (fTS"]) . The 
Eq. (fTS"]) contains the QCD evolution matrix U. Since the 
U matrix is modified due to the removable singularity for 
/ = 3 as given in Eqs. p7H40p . the QED part, the R ma- 
trix, will also change correspondingly. In the following 
we provide this modified version of the R matrix. As we 
mentioned in the previous section, we can not use the S 
and J matrices because they are singular. Instead, we 
need to define a new matrix which is finite: 



Sij(l - 5i,8<>j,7) 



(44) 



Note that the singular part is subtracted away so that 
the H matrix is finite. 



5 



We now express the R matrix as follows: 



where is 



2tt 

i?(mi,m 2 ) = — — VK(mi, mi)V~ 
Po 

K(m 1 ,m 2 ) = K^°\mi,m 2 ) 



1 - 



(45) 



(46) 



where the leading term is given in Eqs. (|28H30p . The 
matrices are given as 



M (2) = y-lp(l)y 



3 



[A' 



4 ]ij 



a s ( mi )[H K(% 



,(m2)[KW H] 



Sis 



2po 



[h]ij +S 3 7 



M iS ^87 



2/?o 



(47) 
(48) 
(49) 

(50) 



(51) 



A complicated calculation lead to the following expres- 
sions for the I\ and I 2 matrices: 



[I: 



[I: 



2\ij 



a s (m2) 
ai—aj \ a s (mi) 



1 / a s (m 2 ) 

2 \ q s (mi ) 



_____ _s ___) 

ai—aj \ ce s (mi) 



In 



1 [ a a (m 2 ) 
' 2 I o.(mi) 



a s (m 2 ) ) (at—aj)' 2 
2 

In [ £4124 

ce s [mi) 



a s (rii2) 



a s (mi) y (a.j —aj ) 2 
fl s (mij 



Q 3 (m 2 ) 
a s (mi ) 



a s (m 2 ) 
a s (mi) 



a a (m 2 ) 
a s (mi ) 



a,(m 2 ) 
a s (mi ) 



if a t ^ aj 
if di = aj 

if a; 7^ aj 
if a,- = cla 



(52) 



(53) 



r 



Note that the leading order contribution is the same 
as before. The only change is localized in the ma- 
trices. In particular, the matrices for (i = 1, 2, 3) 
corresponds to the matrices once we substitute the 
S matrix by the H matrix. The K± matrix represents 
the contribution from the removable singularity of the 5* 
matrix. The key point is that the if} matrices are fi- 
nite and can be used numerically, whereas the are 
divergent. 



QED part of the RG evolution matrix are finite and can 
be used for numerical studies with Nf = 2 + 1. These 
results are currently being used to analyze the data sets 
of the staggered e'/e project (la ]. In fact, this work is a 
part of that project. 



CONCLUSION 
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APPENDIX: ALTERNATIVE DERIVATION OF 
THE PURE QCD EVOLUTION MATRIX AT NLO 

In this appendix we give an alternative derivation of 
the result for the pure QCD evolution matrix U (mi, m^) 
at NLO in the singular case (Eq.'s (J3TJ) and [|4"0])), Since 
the dependence of U (mi, 7712) on m% enters through 
a s (mi) = g(rrii) 2 /4n , we will sometimes denote the evo- 
lution matrix by U(gi,g 2 ) when it is convenient. Follow- 
ing the earlier approach of Buras et al. [!, [I?) , we express 
the evolution matrix as 

U(g,9o) = ( l +^2 J (9)) U m (9,9o) (1+^ J (9o] 

(A-l) 

and seek to determine J(g) from the differential equation 
which characterizes U(g,go): 



„(°) 



7^°' vanishes, 



In the singular case where 2/3q + 7, 
(|A.3|) does not admit a consistent solution at lowest order 
based on an expansion J(g) = J + O(g). This case was 
left untreated in the work of Buras et al. However, this 
case can be readily dealt with after realizing that J{g) 
need not reduce to a constant matrix J for g — > 0. In fact, 
all that is required is that g 2 J(g) — ► for g — > 0, so that 
the expression (IA.1[) for U(g,go) reduces to U^°'(g,go) 
for small g as it must. This allows for the possibility 
that J(g) diverges for g — > 0. Indeed, in the singular 
case (which can only occur when i 7^ j) (|A.3|) reduces at 
lowest order to 



Po g Sy (5) 
which has the solution 



G 

^iCflO = — ^ lo g(5) + Ci 
Po 



(A.5) 



(A.6) 



(A.2) 



Substituting (|A. 1|) into (|A.2|) . and introducing S(g) via 
J( 5 ) = VS(g)V~ 1 as in Eq. (J2U, we find the following 
differential equation for the matrix elements of S(g): 



Po g S' l3 (g) + (2/3 + 7 | 0) - 7 f )S ij (g) 
° {£h { $-G + Otf))S{g)) 

_ ft (0) 



167T' 



where is an undetermined integration constant. In 
fact we do not need to determine cy since it turns out 
not to contribute to the evolution matrix at NLO. To 
see this, recall that the evolution matrix is determined 
at NLO by A(mi,m2) as in Eq. (|37|) , where now 

A(mi,m 2 ) = a s (mi 

-a s (m 2 )V- 1 U( Q) (m 1 ,m 2 )VS{m 2 ) . 

(A.7) 



Po 



ll Sij ~ Gij + 0(g ) (A.3) Re-expressing (TOt as 



To solve this equation at leading order in g it is necessary 
to consider separately the cases where 2po + 7^ — 7] 
is non- vanishing and vanishing. In the former case, a 
consistent solution is obtained at lowest order by taking 
J(g) = J + 0(g) where J is a constant matrix. Then 
S(g) = S + O(g), and the lowest order part of (|A.3|) 
becomes 



Sij(m) 



G 

-^log(a s (m)) + c' y - 



(A. 



(2/?o + 7) 



(0) 



Dividing by 2/3 + 7}° — 7]°^ gives the expression for SV, 
stated in Eq. ([22L which was the one obtained previously 
by Buras et aL[£, [T3| ■ 



Po 



Gi 



(A.4) 



(where c[j — Cij — jj^ log(47r)) and substituting this into 
the expression for A(mi,m 2 )ij obtained from (|A.7|) (re- 
calling from Eq.'s ([15 ]) - pO] ) that y _1 [/°(mi, m 2 )F is di- 
agonal) we easily find that the constant c^- drops out 
and our previous expression Eq. (|40p is reproduced. This 
completes the alternative derivation of the NLO expres- 
sion for U(mi,m 2 ). The argument also shows that for 
calculations involving U(mi,m 2 ) at NLO we may take 
Sij in the singular case to be given by (|A.8[) wit h cj. - = 0. 
This is useful for deriving the expressions Eq.'s (|50[)-f 
for the QED part of the evolution matrix at NLO. 
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